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Abstract — The properties of backward-wave second harmonic meta-reflector operating in pulse 
regime are investigated. We show that ratio of the income pulse length to the thickness of 
the metaslab determines its important operational properties which are in contrast with those 
manifested at common second harmonic generation settings. 


1. INTRODUCTION 

Light can be described as travelling electromagnetic wave (EMW) with phase velocity v p h directed 
along the wave vector k whereas its energy flux is represented by the Poynting vector S: 

Vp/i = (k /k)(c/n), S = (c/47t)[E x H] = (c 2 k/47rwe)iL 2 = (c 2 k/47ru;/r)E 2 . (1) 

Here, c is speed of light, n is refractive index of the propagation medium, E and H are electric 
and magnetic components of the wave. Refractive index is determined by electric permittivity e 
and magnetic permeability n of the medium at corresponding frequency. In all naturally occurring 
materials, e > 0, /j > 0, n = y/jli and vectors S and k are co-directed. In the case of e < 0 
and n < 0, refractive index becomes negative, n = — y/JIe, and vectors S and k - contra-directed. 
Such waves are referred to as backward EMW (BEMW). The possibility to produce backward light 
became achievable only recently owing to the advent of nanotechnology and optical metamaterials 
which has led to revolutionary breakthrough in the concept and in the numerous applications of 
the linear optics [T|. Extraordinary coherent nonlinear optical (NLO) processes were predicted 
in for the cases of coupled ordinary and BEMW waves. Among them are second harmonic 
generation (SHG) pH!], optical parametric amplification and frequency-shifted nonlinear reflectiv¬ 
ity @H6]. Metamaterials (MM) are artificially designed and engineered materials, which can have 
properties unattainable in nature. Current mainstream in fabricating negative-index MM (NIM) 
slabs relies on engineering of LC nanocircuits - plasmonic mesoatoms at the nanoscale with negative 
electromagnetic response. Extraordinary coherent nonlinear optical frequency-converting propaga¬ 
tion processes predicted in NIMs have been experimentally realized to date in the microwave [7] 
and in the multilayered plasmonic optical MM [8]. However, a different, more general approach, 
to engineering MM that can support BEMW is possible Pino]. It is based on the fact that in a 
loss-free isotropic medium, energy flux S is directed along the group velocity v g : 

S = v g U, v g = grad k w(k). (2) 

Here, U is energy density attributed to EMW. It is seen that group velocity is directed against 
the wavevector if dispersion dui/dk becomes negative. Basically, negative dispersion dco/dk < 0 
can appear even in fully dielectric materials with particular composition of its structural elements. 
This opens an entirely novel research and application avenue. Various particular realizations of 
negative dispersion were proposed so far HM5]. Basically, many hyperbolic MM and specially 
designed waveguides can support BW electromagnetic modes, see, e.g., [T6HI8] . This opens new 
avenues for creation of BW photonic devices with unparalleled functional properties. As regards 
coherent nonlinear-optical propagation processes, critically important is to provide for coexistence 
of coupled ordinary and backward waves which frequencies and wave vectors satisfy to energy and 
momentum conservation law (phase matching). Such possibilities were described in ref. [131415] . 

This paper is to investigate unusual dependence of BW SHG in the pulse regime. We show that 
conversion efficiency and properties of generated SH pulses depend on the ratio of the input FH 
pulse length to thickness of the NLO MM slab. Losses are included in the consideration. 


2. BASIC EQUATIONS 

As noted, properties of SHG will be investigated for the cases where one of the coupled waves is 
ordinary and the other one is BW. To achieve phase matching, wave vectors of the fundamental and 
the SH waves must be co-directed. This means that in BW setting pulse of the SH will propagate 
against the pulse of fundamental radiation. Corresponding basic equations are as follows. Electric 
and magnetic components of the waves and corresponding nonlinear polarizations at uq and at 
oj 2 = 2uq are defined as 


{£, H}j = Re {E , H}j exp {i(kjZ — ujjt)}, {V, M}j = Re {P, M} J exp {ifkjz — ujjt)}, 

{P,M}? L = E,H}\,{P,M}$ L = xil 2 {E,i7}?, 2 X g^ = 

Amplitude E\ of the first harmonic (FH) and of the SH, E 2 , are given by the equations: 
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Here, Vi > 0 and ai j2 are group velocities and absorption indices at the corresponding frequencies, 
Xeff = Xe 2 effective nonlinear susceptibility, A k = k 2 — 2k\. Parameter Sj = 1 for ordinary, and 
Sj = —1 for backward wave. With account for k 2 = n 2 (u;/c) 2 , n\ = n 2 = s 2y Je 2 [i 2 , we 

introduce amplitudes e 3 = y/\ej\/kjEj, aj = ei/eio, coupling parameters se = \/k 1 k 2 /l€ 1 e 2 l47rx ( e ff 

and g = aeEio, loss and phase mismatch parameters 5 i j2 = ai )2 L and A k = A kl, slub thickness 
d = L/l, position £ = z/l and time instant r = t/ At. It is assumed that Ej$ = Ej(z = 0), 
l = v\A t is the pump pulse length, Ar is duration of the input fundamental pulse. Quantities |aj| 2 
are proportional to the time dependent photon fluxes. Then Eqs. m and d6|) are written as 
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In the case of magnetic nonlinearity, 
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equations for amplitudes Hj take the form: 
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Then, equations for amplitudes aj = rrii/mio, where rrij = y/\fij\/kjHj, take the form of ((TJ) and 
(f8l). with coupling parameters given by ae = \/k\k 2 /1 g 1 n 2 \and g = seH w . In both cases, 
Zo = g _1 is characteristic medium length required for significant NLO energy conversion. Only 
most favorable case of the exact phase matching A k = 0 will be considered below. 


3. EXTRAORDINARY PROPERTIES OF THE BACKWARD-WAVE SECOND 
HARMONIC GENERATION IN THE PULSED REGIME 

Major difference in SHG in ordinary, positive index materials (PIM), and BW (NIM) materials 
in the continuous wave (cw) regime is summarized and illustrated in Fig. |T] [3ll4l ll9| . In order 
to achieve phase matching, wave vectors for FH and SH must be co-directed, which dictates co¬ 
directed energy fluxes in PIM [Fig. Oja)] and contra-directed in NIM [Fig. [l](b)] . In a PIM, photon 
flux in FH depleats whereas photon flux in SH grows across the medium so that the latter one 
may exceed that in FH by the exit of the NLO slab. In lossless PIM, sum of the photons in SH 
and of the photon pairs in FH is conserved in any point of the NLO media. On the contrary, FH 
and SH propagate in the opposite directions in a NIM. Consequently, the difference of the above 
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Figure 1: Fundamental difference in the properties of continuous-wave phase matched SHG in a loss-free 
ordinary and BW materials, (a) and (b): Difference in the NLO coupling and photon fluxes geometries, 
(c) and (d): Photon fluxes reduced by the input magnitudes across the slab. Here, the descending dash- 
dotted blue line depicts FH flux Si and the ascending dotted blue line depicts S 2 for SH, both correspond 
to ordinary waves in a positive-index (PIM) material. The descending solid red line depicts Si flux for 
backward-wave FH, whereas the descending dashed red line depicts ordinary-wave S 2 flux for SH, both in 
a frequency double-domain negative-index (NIM) slab, (c) gL = 0.5, (d) gL = 3. (e): Output transmitted 
fundamental (the descending blue dash-dotted line) and SH (the ascending blue dotted line) fluxes at z = L, 
both are ordinary waves in a PIM. The descending solid red line displays transmitted BW flux Si at z = L 
and the ascending dashed red line depicts output ordinary SH flux at z = 0, both correspond to a NIM. 

indicated numbers is conserved so that the number of photons huj 2 at 2 = L is always equal to 
zero and is less than the number of photon pairs in FH at z = 0 [Fig. [D(c,d)]. The rate of the 
changes across the slab depends on the field strength in the input FH beam [cf. Fig. [ljc) and (d)]. 
Figure [0(e) depicts dependence of the photon fluxes in the FH at the exit of the NLO slab and for 
SH at the exit, z = L, for the PIM slab and photon flux of SH from the opposite , z = 0, edge of 
the NIM slab. Obviously, the outlined properties of SHG in the case of ordinary FH and BW SH 
are similar. Figures [D(a-e) exhibit fundamental differences in the the properties of the SHG which 
involves backward light as compared with the SHG in ordinary crystals. 

Unparalleled properties of BWSHG in the pulsed regime stem from the fact that it occurs only 
inside the traveling pulse of fundamental radiation. Generation begins on the leading edge of the 
pulse, grows towards its trailing edge, and then exits the pulse with no further changes. Since the 
fundamental pulse propagates across the slab, the duration of the SH pulse is expected to be longer 
than that of the fundamental one. Depletion rate of the FH radiation across its pulse length and 
the conversion efficiency must depend on its initial maximum intensity, phase and group velocity 
matching. Ultimately, the overall properties of BWSHG such as the output pulse length and the 
photon conversion efficiency can be foreseen dependent on the ratio of the fundamental pulse and 
slab lengths. Investigation of the indicated dependence is the major goal of this work. 

The input pulse shape is chosen close to a rectangular form 

F(t) = 0.5 (tanh T ° ^ —— — tanh ^ , (11) 

V or or J 


where 5t is the duration of the pulse front and tail, and tq is the shift of the front relative to t = 0 . 
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Figure 2: Backward-wave second harmonic generation in the pulse regime. The red dash-dotted plot depicts 
input (at 2 = 0) pulse of FH and the red dashed plot depicts transmitted pulse at z = L. The blue solid 
plot depicts output contra-propagating pulse at doubled frequency at 2 = 0. r = t/ At, gL = 3, d = L/l. 
(a, b): losses are neglected, (c, d): u\L = 2, a 2 L = 0.2. (a,c): d = 0.1. (b,d): d = 0.5. S'i(O) = 0.9900. (a) 
S'i(L) = 0.0934, S 2 (0) = 0.4485. (b) Sx{L) = 0.1117, S 2 (0) = 0.4391. (c) S^L) = 0.03, S 2 { 0) = 0.2865. (d) 
Si(L) = 0.0343, 5 2 (0) = 0.2695. 


Parameters 5t = 0.01 and to = 0.1 were selected for numerical simulations based on the system of 
partial differential equations (J7]) and (0). The results are illustrated in Figs. [2] and [3l A rectangular 
shape of the fundamental pulse, Tf = |ai(2)| 2 /|aio| 2 , is shown at 2 = 0 when its leading front enters 
the medium and the results of numerical simulations for the output fundamental pulse when its 
tail reaches the slab boundary at z = L. Shape and conversion efficiency of the output SH pulse, 
772 = |a 2 ( 2 )| 2 /|aio| 2 , traveling against the z-axis is shown when its tail passes the slab’s edge at 
2 = 0. Here, phase and group velocities of the fundamental and SH pulses assumed equal. Figure 
Ha) corresponds to the fundamental pulse of ten times and Fig. 0L) of two times longer than the 
slab thickness. Losses are neglected. It is seen that pulse shapes of SH are significantly different 
in (a) and (b) cases. It is because the transient periods when only part of the fundamental pulse 
is inside the slab is longer relative its duration in the second case. Since the metaslab is lossless, 
the outlined properties satisfy to the conservation law: the number of annihilated pair of photons 
of in the FH (S 10 — Sil)/ 2 per pulse is equal to the number of output SH photons S 20 ■ Here, 
Si, 2 Z are numbers of photons integrated over the corresponding pulse duration and normalized 
by the corresponding integrated input numbers. Note, that maximum possible photon conversion 
efficiency for SH is 0.5. Figures H c id) show effect of losses which roughly correspond to one of the 
model of the metaslab made of free-standing carbon nanotubes P3HI5]- 

































Figure 3: Input, transmitted and reflected pulses. (a,c): d = 5, (b,d): d = 10. <Si(0) = 0.9900. (a) 
Si(L) = 0.3607, S 2 ( 0) = 0.3111. (b) Si(L) = 0.5116, £ 2 ( 0 ) = 0.2403. (c) Si(L) = 0.0734, S 2 ( 0) = 0.1273. 
(d) Si(L) = 0.0923, 52(0) = 0.0787. Notations and other parameters are the same as in Fig. [2] 


Figures [3] (a,b) illustrate the effect of further pulse shortening. Here, Fig. [3] (a) displays input 
pulse five times and Fig. [3] (b) - ten times shorter than the metaslab thickness. From comparing 
figures [2] and [3] one concludes that difference in the pulse lengths grows so that SH pulse becomes 
much longer compared to the FH pulse. Counter-intuitively, quantum conversion efficiency per 
pulse decreases with decrease of the input pulse duration even though its maximum field strength 
remains unchanged and losses change SH pulse shape. 

4. CONCLUSIONS 

Properties of frequency doubling nonlinear-optical metamirror operating in the pulsed regime are 
investigated. It is proposed to be engineered based on the metamaterial which supports backward 
light waves at one of the frequency and ordinary wave at another frequency. It is supposed that 
phase velocities can be adjusted equal and are co-directed, whereas energy fluxes of the incident 
and generated second harmonic light are contra-directed. References are given to the work that 
prove such a possibility. Physical principles underlying fundamental difference in the properties 
of second harmonic generation in the proposed and standard settings are discussed. A set of 
partial differential equations which describe such a reflector with the account for losses are written 
and solved numerically. It is shown that unlike second harmonic generation in standard settings, 
contra-propagating generated second harmonic pulse may become much longer then the incident 
fundamental one and the difference grows with decrease of the input pulse length as compared with 


























thickness of the metaslab. The revealed properties may manifest themselves beyond the optical 
wavelength range. 
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